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ABSTRACT 

Gravitational lensing observations of massive X-ray clusters imply a steep character- 
istic density profile marked by a central concentration of dark matter. The observed 
mass fraction within a projected radius of 150 kpc is twice that found in state-of- 
the-art dark matter simulations of the standard ACDM cosmology. A central baryon 
enhancement that could explain this discrepancy is not observed, leaving a major 
puzzle. We propose a solution based on the merger histories of clusters. A significant 
fraction of the final dark matter content of a cluster halo originates within galaxy- 
sized halos, in which gas can cool and compress the dark matter core to high densities. 
The subsequent tidal stripping of this compressed dark matter occurs in denser re- 
gions that are closer to the center of the cluster halo. Eventually, the originally cooled 
gas must be dispersed into the intraclustcr medium through feedback, for consistency 
with observations that do not find central baryon enhancements in clusters. Still, the 
early adiabatic compression of the galactic dark matter leaves a net effect on the clus- 
ter. Using a simple model for this process, we show that the central cluster profile is 
substantially modified, potentially explaining the observed discrepancy. 
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1 INTRODUCTION 

Recent observations have confirmed our basic understanding 
of cosmology and showed an impressive consiste ncy with the 
predi c tions of the standard ACDM model (e . g..lAstier et al.l 
l200d : iPercival et "ail 120071 ; iKomatsu etatl 120091 ). In this 
model, a cosmological constant dominates the cosmic mass 
budget today, but galaxies and other structures were assem- 
bled earlier, primarily out of cold dark matter. This medium 
of non-interacting, low velocity-dispersion particles, started 
out with small Gaussian density perturbations that were 
subsequently enhanced by gravity. While the model success- 
fully matches observations of the large-scale anisotropies of 
the cosmic microwave background and the large-scale struc- 
ture in galaxy surveys, it is also important to test its validity 
on smaller scales. The abundance and structure of non- linear 
objects are potentially sensitive probes of the properties of 
dark matter (e.g., whether it is cold) and of the density 
fluctuations (e.g., whether they are Gaussian). However, gas 
cooling and astrophysical feedback complicate the interpre- 
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tation of observations regarding the dark matter distribution 
in galaxies. Thus, it is most attractive to study the largest 
virialized objects, namely X-ray clusters, in which most of 
the gas is too hot and rarefied to cool and is thus expected 
to trace the gravitational potential. 



The mass profiles of galaxy clusters can be measured 
directly through gravitational lensing. Observations of the 
most massive clusters now find dozens of multiply-imaged 
background sources, allowing a precise measurement of the 
central 2-D mass distribution in each cluster as projected 
on the sky. Also crucial for characterizing each cluster is its 
total virial mass, which can be measured precisely by supple- 
menting the central strong lensing signal with weak lensing 
distortions measured out to the cluster edge. It is useful to 
characterize the total, projected profile with one scale, the 
effective Einstein radius te (or angle #e) defined so that a 
circle of that radius around the cluster center contains a 
mean enclosed surface mass density E equal to the critical 
density for lensing, E cr = [c 2 /(47tG)]Z>os/(-Dol-Dls), where 
D denotes various angular diameter distances (Observer- 
Source, Observer-Lens, and Lens-Source). This definition is 
motivated by the Einstein ring radius of an axisymmetric 
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lens, but is nonetheless a useful measure of the central mat- 
ter content even for asymmetric clusters. 

Current N-body simulations of galaxy clusters in 
ACDM produce samples of thousands of halos with virial 
mass M v i r > 10 14 A/q. These simulations are becoming suf- 
ficiently large and detailed to yield the predicted spread in 
cluster halo parameters, and to allow a quantitative assess- 
ment of the inherent bias in observing cluste rs in projection 
and selecting them by lensing cross-section l|Hennawi et al.l 
120071 : iNeto et all 120071 ). In general, the density profiles of 
the simulated clusters are relatively shallow and seemingly 
at odds with recent c areful lensing studies of massive clusters 
dKneib et al.l 120031; ICayazzi et all 12003); iBroadhurst et al ' 



model. Subse quent work has ge n erally supported this con- 
clusion (e.g.. loguri et al.| 120091 : lOguri fc Blandfordl |2009l ; 



20051; iKlinget al. 2005; Limousin et al.ll2007l; iBradac et al 



2008; Halkola et al. 120081 ; Umetsu fc Bro adhurst 2008). 

Th is discrepancy was recently hi ghlighted and quanti- 
fied by IBroadhurst fc Barkanal (|2008l ) , who carefully com- 
pared observations of four well-constrained massive clus- 
ters to the predictions of the numerical simulations. They 
emphasized the importance of comparing directly the pro- 
jected 2-D mass distributions in the observations and the 
simulations, using the virial mass (M v ir) and the effective 
Einstein radius (ve) as two numbers that characterize the 
degree of concentration in each cluster halo. They derived 
the theoretical predictions for cluster lensing in ACDM by 
sta rting with t h e dist ribution of 3-D halo profiles measured 
by INeto et all (|2007h in the Millennium simulation, and 
th en correcting it for le nsing and projection biases based 
on lHennawi et all l|2007h . Comparing the resulting distribu- 
tion with the observed te for four clusters - A1689, C10024, 
A1703, and RXJ1347 - and including the expected spread 
in profiles as well as the measurement errors, they found 
that each cluster was discrepant at the 2— a level (all with 
an unusually large te giv e n M v j r ), yielding a combined 4-cr 
discrepancy. iDuffv et al.l l|2008h recently found that simu- 
lated cluster concentrations are even lower when using the 
most updated cosmological parameters (which h ave a lower 
power spectrum normalization than assumed bv lNeto et al.l 
((2003)), though the effect for the most massive clusters is 
only at the level of ~ 10%. 

IBroadhurst fc Barkanal (|2008T ) suggested that gas 
physics is unlikely to affect significantly the Einstein radius 
of a cluster. This radius of ~ 150 kpc is typically observed to 
enclose a projected mass of ~ 2xl0 14 M Q , or a mass of ~ 1 x 
1O 14 M0 within the same 3 -D radius. Using the sim ple model 
of adiabatic compression (|Blumenthal et al]|l986l ). they es- 
timated that the observed 3-D mass within the Einstein ra- 
dius can be obtained if gas cooling increases the enclosed 
baryonic fraction within this radius to ~ 1/3, twice the cos- 
mic baryon fraction. Indeed, hydrodynamic simulations pro- 
duce clusters that are as centrally concentrated as those ob- 
served, likely due to their "overcooling" problem which pro- 
duces just such an increa se in the central baryon fraction, 
with most of it in stars (|Kravtsov et al.l 120051 ; iRozo et al.l 
2008). An increase of this sort apparently does not occur 
in real clusters, where the baryonic (gas+stellar) mass frac- 
tion within the Einst ei n radius is below the cosmic value 
Lin et al.l |2004|; iLaRoque et al.l |2006| ; iLemze et all 
' Afshordi et al.l 120071: IVikhlinin et all 120091 '). Thus. 
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Broadhurst fc Barkanal (|2008f ) concluded that cluster halo 
profiles present perhaps the clearest, most robust, cur- 
rent conflict between observations and the standard ACDM 



IZitrin fc Broadhurst! I200E ) , though only a large unbiased 
cluster sample with precise strong and weak lensing mea- 
surements would be completely conclusive. 

In this paper we propose a novel process that could re- 
solve the apparent discrepancy between cluster observations 
and existing ACDM simulations. A significant fraction of the 
final dark matter content of a cluster halo originates within 
galaxy-sized halos, in which gas can cool and compress the 
dark matter core to high densities. In § 2 we develop a simple 
model for this adiabatic compression and for the subsequent 
tidal stripping of the dark matter within the cluster halo. 
We then show in § 3 that the central mass profile of the 
galaxy cluster is substantially modified by compression of 
the galactic halos that it swallows, even if the cooled galac- 
tic baryons are later redistributed throughout the cluster. 
We show quantitatively that this can potentially explain the 
observed discrepancy. Finally, we summarize our conclusions 

and caveats in § 4. 

We note that lMaccio et all {2006) found in a simulation 
of a Galactic halo that cooling produced a more highly con- 
centrated dark matter profile for the host halo, but there, as 
in the cluster simulations, this may have been due directly to 
adiabatic compression in response to the large concen t ration 
of baryons in the center of the host halo. iDolag et al.l (|2009h 
included radiative cooling and stellar feedback in simulations 
of cluster formation, and while they focused on the structure 
of the galactic subhalos, they also found a ~ 15% increase in 
the effective concentration of the cluster halo; this, however, 
may still be partly due to an increased baryon content near 
the cluster center. Also. ISommer-Larsen fc Limousin! (|2009h 
artificially eliminated the over-cooling problem in their clus- 
ter simulations and found only a very minor effect of the 
baryons on the total mass profile in this case; however, their 
artificial scheme may have also eliminated the effect we an- 
alyze. Our analytical approach allows us to cleanly separate 
the effect of density-enhanced tidal stripping from a simple 
overall adiabatic compression of the cluster halo. 



2 MODEL 



We a ssume the standard ACDM cosmology ()Komatsu et al.l 
2009), with a dimensionless Hubble parameter h = 0.7 
and density parameters Q m ~ 0.28 (dark matter plus 
baryons), = 0.72 (cosmological constant), and fit, = 
0.046 (baryons). We also denote the cosmic baryon fraction 
by f b = (Q b /n m ) = 0.16. 

Consider a halo that virialized at redshift 2 in a flat 
ACDM universe. The critical density at z is 



= SE[<Ui + z) 3 + «a] 



(1) 



Numerical simulations of hierarchical halo formation indi- 
cate a roughly univ ersal spherically-aver aged density profile 
for virialized halos (|Navarro et al J ll997l , hereafter NFW): 

P( r ) = Pc 7- 1 72 ' ( 2 ) 



where the radius r is divided by the scale radius r s 



Tidal Enhancement of Cluster Concentrations 3 



TVir/cvir with r v i r being the virial radius, and the charac- 
teristic density 5 C is related to the concentration parameter 
Cvir by 



Sr.= 



vii 



3 ln(l + c vir ) - c vir /(l + c vir ) ' (3) 

where A c is the virial density in units of p z c . For a halo of 
virial mass A/ V j r at a given redshift z, the profile is fully 
specified by the parameters A c and c v ir. We adopt the con- 
vention of a fixed A c = 200 at all redshifts, for consistency 
with the simulation analyses whose results we use. 

We adopt a simple analytical model of tidal stripping 
that, in particular, has been previously used to understand 
how the density profiles of satellite sub-hal os produce the 
centr al NFW profile of their final host halo (|Sver fc White! 
1998). In this model, material originally at a radius £ within 
the satellite ends up, after tidal stripping, at a radius r 
within the host halo so that the mean enclosed satellite den- 
sity within £ equals the mean enclosed host density (before 
the stripping) within r: 



PhostM = Psat(£) 



(4) 



This condition corresponds to a resonance in dynamical fre- 
quencies between the circular orbit of a mass element at 
£ around the satellite and that of a circular satellite's or- 
bit around the host halo, inducing an energy transfer that 
strips the mass element from the satellite. Setting the host 
tidal force equal to the internal gravitational force within 
the satellite yields a very similar stripping radius. 

We combine this strip ping model with the sim ple model 
of adiabatic compression (Blu menthal et alj fl986) in which 
conservation of angular momentum implies that the quantity 
rM(r) (assuming spherical symmetry) is fixed. We assume 
that both the host and each stripped satellite start out with 
NFW profiles. In satellites within the mass range of galax- 
ies, the baryons cool and condense to the center, inducing 
a change in the surrounding dark matter halo. Specifically, 
adiabatic compression moves a mass shell initially at , con- 
taining a mass M aa t(£i), to a final radius £f = £iM 8a t(£x) / 'Mf , 
where the final enclosed mass is larger by a factor 



1-fb + fb 



M si 



(•») 



M sat (&) J " J "M s ^i) ' 

where M sat without an argument denotes the total virial 
mass of the satellite. The stripping model then implies that 
the same mass shell ends up at a radius r{ in the host cluster 
halo, where 



Phost(rf) = Psat(Ci) X 



Mi 



M sat (£i 



(6) 



in terms of the initial enclosed density p S at(£i) in the satellite 
(i.e., before the adiabatic compression). The power of four 
on the right-hand side results from the increased mass (one 
power) and decreased radius (hence three powers in the den- 
sity). Note that the cooled baryonic cores of the galaxies are 
much denser than their surrounding dark matter halos and 
so we have assumed that the baryonic cores are not tidally 
stripped. 

In deriving equation ([5]) we have assumed that in galac- 
tic halos, where the cooling time of the virialized gas is much 
shorter than the Hubble time, the full baryonic content of 
the halos initially cools toward the center of the halo, con- 



densing the surrounding dark matter before it gets stripped. 
Observationally, even the fraction of galactic halo baryons 
that are in stars today is not well known, since the hot 
baryons in halos are difficult to detect, while total masses of 
galactic halos are difficult to measure accurately (and can 
be used to obtain the total gas mass only with the added 
assumption of a halo baryon fraction that equals the cosmic 
mean) . The best estimate for the total mass of stars and stel- 
lar remnants today, as a fraction of the total baryonic mass 
that lies within virialized reg ions of galaxies, is ~ 10%. with 
an uncertainty of order 50% dFukugita fe Peeblesll2004l ). For 
our own Milky Way galaxy, the disk and bulge may mak e 
up as much as 40% of the halo baryons (|Xue et alj|2008l ). 
Regardless of the precise fraction today, it is plausible to 
assume that most of the gas in galaxies initially cooled and 
was later expelled over time from the central region back 
into the halo, through supernova or quasar feedback. Thus, 
we expect that if galactic halos were stripped within the 
cluster relatively early, then star formation and feedback 
did not have much time to operate prior to the stripping. 
This scenario does not conflict with the fact that clusters 
only virialized relatively recently, since we are focusing here 
on the stripping that formed the inner regions of clusters, 
within a tenth of the virial radius, and this likely occurred 
long before the entire cluster virialized. 

In clusters, the originally cooled gas must eventually 
have gotten dispersed into the intracluster medium, since 
observations find only a small fraction of cluster baryons 
residing near the center. In the cluster environment, this 
gas redistribution can be facilitated by interactions among 
galaxies or with the intracluster medium, in addition to in- 
ternal galactic feedback. For simplicity we assume that the 
final baryon distribution is similar to that of the dark mat- 
ter, i.e., the final baryon fraction is uniform and equal to the 
cosmic value. Within our model above, this effectively means 
that the satellites contribute only M sa t(fi) to the mass en- 
closed within the final cluster radius rt (and not the full 
Mi). 

In this picture, when the cooled galactic baryons get 
redistributed throughout the final cluster halo, they may 
cause a partial adiabatic expansion of the halo. However, this 
should roughly cancel the initial adiabatic contraction of the 
surrounding halo material when the satellites enter and are 
stripped; thus, we neglect both the initial halo contraction 
and the later expansion, as we do not expect a significant 
net effect. This is different from the main process that we 
focus on, where the early adiabatic compression causes the 
satellites' dark matter to be stripped at smaller cluster radii 
than it would otherwise, leaving a net effect on the cluster 
in the end. 

The stripping model that led to equation (|5J assumes 
that the host halo is dominant, and that the satellites con- 
tribute only a small fraction of the enclosed mass at the 
stripping radius. This assumption breaks down, however, 
when the satellites dominate, and we do expect such a 
regime; indeed, it is plausible (and indicated by pure dark 
matter simulations) that the dense core of the cluster halo 
arises entirely from the original cores of the accreted satel- 
lites, since only l ower-density material i s accreted later onto 
the cluster halo (|Loeb fc Peebles! 120031 ; iGao et al.ll2004l ). In 
this high-density regime where the satellites dominate rela- 
tive to the pre-existing halo material, we expect the satellite 
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cores to simply settle in the host core, preserving in the end 
their original densities, since they first adiabatically com- 
press and later re-expand. This limit is consistent with the 
equal-density relation in equation Q; while this equation 
was originally derived for stripping onto a dominant host 
halo, we can use it to capture the dominant satellite limit, 
if we obtain from it the radius n at which the satellite mass 
ended up (reinterpreted as contributing part of the halo pro- 
file rather than adding mass on top of p\\ ost )- We leave for 
future work a detailed analysis of the complex transition re- 
gion between the two limits of equation Q (density preser- 
vation) and equation (density enhancement), and here 
we adopt a simple interpolation between them. This is rea- 
sonable given our limited goal of examining whether the 
halo profile can be substantially modified at all. Thus, we 
determine the stripping radius by solving (for a given £i): 



Phost(n) = p S at(6)x < / 8a t(r f ) + [1 - / 8a t(r f )] 



Ms 



M sat (£i) 



,(7) 



where / 8a t(ff) is the fractional contribution of the satellites 
to the enclosed mass within the stripping radius rf, and 
the right-hand side of this equation interpolates between 
equation (valid in the limit / sa t(ff) = 1) and equation ((6]) 
(valid when / sat (rf) = 0). 

In order to solve the model, we must determine various 
masses. First, based on the host NFW profile in the absence 
of cooling and adiabatic compression, we obtain the enclosed 
host mass at each radius, M^ ost (rf) = |vrrf Ph ost (rf). We use 
equation @ to calculate the fractional contribution of the 
satellites to this enclosed mass in the absence of cooling, 
/sat^f). Adiabatic compression and stripping then replaces 
the satellite contribution by Mf, and thus in equation (J7J) 
we set 



/sat(r f ) 



(8) 



After solving for rf, the final enclosed mass within this ra- 
dius, after feedback redistributes the baryons, is 



M bost (r f ) = [1 - / s ° at (r f )]M h ° ost (r f ) + M sat (£) 



(9) 



Within our simple model, the redistribution of mass 
depends only on the density profile of the satellites, i.e., on 
the functional form of M aa ,t (£i)/M sa .t versus p sa t(£i), and not 
on the number of satellites or their individual total masses. 
Since we are interested in the galactic sub-halos that end 
up in the cluster, we can adopt the typical value of accre- 
tion redshift z sa t and NFW concentration c v j r for such halos, 
and effectively calculate stripping of one satellite that con- 
tributes some total fraction / sat of the host cluster mass. 
Within the model, this single satellite represents the cu- 
mulative effect of all the individual galactic satellites that 
merged into the final cluster halo. 

In order to quantify the effect of compressed galactic 
halos, we must estimate the effective / sa t, i.e., the fraction 
of cluster dark matter that arrived from within galactic ha- 
los in which the baryons were able to cool. We can obtain a 
theoretical estimate for the fraction of cluster dark matter 
that passed in its merger history through galactic halos, on 
the way to becoming part of the final cluster halo. For con- 
creteness, let us consider the progenitor distribution at var- 
ious redshifts of a 10 15 M« cluster halo at z ~ 0.2 - 0.4. The 



extended Press-Schechter model l|Bond et al.ll 199lh then im- 
plies that ~ 25% of the cluster mass was in halos with masses 
in the range 10 10 - 10 12 Af Q at z ~ 2.5 (the redshift that 
maximizes this fraction). A full merger tree would give on 
average at least this value since additional cluster mass that 
was outside this halo mass range at z — 2.5 may have passed 
through galactic halos at other redshifts. 

We can also estimate the fraction / sat from observa- 
tions. A first attempt might proceed as follows. Stars make 
up ~ 1% of the total v irialized mass of massive clusters 
(e.g.. lLagana et al] [2008). which correponds to ~ 5% of the 
total baryonic mass. To find the baryon fraction that was 
associated with the galaxies in which these stars formed, 
we must divide by their average star formation efficiency. 
As noted above, the total stellar mass today is ~ 10% 
of the baryon mass with in virialized regions of galaxies 
jFukugita fc Peebles! 12003 ). This suggests that ~ 50% of 
the cluster gas was processed through galaxies, and thus 
also a similar fraction of the cluster's dark matter was con- 
tributed by stripped satellite galaxies (assuming that the 
total baryon fraction of both these galaxies and the final 
cluster is equal to the cosmic mean fraction). However, this 
estimate depends on the uncertain value of the star forma- 
tion efficiency. We can use metallicity measurements to ob- 
tain a more direct estimate of / aa t. The typical metallicity 
of the intracluster medium in cluster s at redshift z ~ 0.3 is 
0.3 - 0.4 of the solar abundance (e.g-. lMaughan et alJ feOQS). 
while massive galaxies (in clust ers or the field) typi cally have 
a solar abundance or less (e.g.. lEllison et alj|2009h . indicat- 
ing that at least 30 — 40% of the cluster gas must have been 
processed in galaxies in order for the ICM to reach its high 
metallicity value. Within our model, higher / sat values lead 
to more highly concentrated cluster halos (see the next sec- 
tion). Taking into account these various considerations, we 
consider / sat valu es in the rang e 20 — 40%. 

We note that lGnedin et all l|2004h showed that the clas- 
sic adiabatic contraction model that we use tends to overes- 
timate the effect of a central baryon concentration on dark 
matter, compared to simulated profiles. The overestimate in 
the mass profile, however, is under 10% at the radii that we 
focus on (~ 0.1r v i r ), thus justifying our use of the simple 
model. The overestimate does increase at smaller radii and 
is ~ 50% at O.Olrvir, implying that the profiles we find below 
are less reliable in the innermost region. 



3 RESULTS 

In this section we quantify the effect that adiabatic compres- 
sion in galactic satellites can have on the final density profile 
of the host cluster halo . We compare our result s to th e four 
clusters considered by iBroadhurst fc Barkanal (|2008), and 
also make use of their results for the theoretical predictions. 
In particular, we adopt the NFW parameters measured by 
iNeto et"aT. (2007) fo r simulated halos, after correcting them 
based on Hennawi et al.l (|2007h to obtain the effective pa- 
rameters for the population of lensing clusters, observed in 
projection; howe ver, we reduce c v j r b y 10% according to the 
recent results of iDuffv et al.l (|2008l ). This yields a median 
c v i r = 5.5 for the most massive clusters, with a 1—a range 
(for the effective projected c v ir) of 4-7.5 (approximately in 
a lognormal distribution). Studies based on large numerical 
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simulations l|Zhao et alj|2003l ; iGao et al J l2008h have found 
for massive halos a relatively weak decline of cw with in- 
creasing redshift, but a more significant decline for galactic 
mass halos, with c v i r ~ 4 for ~ 1O 11 -1O 12 M0 halos at z ~ 2. 
In our quantitative results, we adopt c v ir = 4 and z sa t = 2 
for the satellites, and consider hosts at z = 0.3 with various 
concentration parameters. As noted above, we assume that 
the dense core of the host arises entirely from the satellite 
cores, so for each host c v ir this normalizes the total satel- 
lite fraction / sat of the cluster's virial mass. In particular, 
we consider / sat = 20% (which implies a host c vlI = 6.8), 
/ S at = 30% (c vir = 7.5), or / sat = 40% (c vir = 8.0). These 
host concentration parameters are somewhat high but still 
near or within the expected 1-a range for cluster halos, as 
noted above. 

Figure [T] shows the effect of adiabatic compression on 
the 3-D mass profiles of clusters. The effect is largest at the 
innermost radii, where the satellites contribute a substan- 
tial fraction of the halo mass. This is true even though our 
model in equation Q suppresses the density enhancement 
in this /sat(ff) — * 1 limit; since the enclosed mass changes 
rapidly with radius in the core, even a slight shift in the 
stripping radius has a large effect on the mass profile. For 
the cases considered, where the fractional mass contribution 
by satellites / sat is 20, 30 or 40%, we find that the cluster 
profile at r < 0.1r v i r is substantially modified by adiabatic 
compression in the satellites, even though the baryonic mass 
is assumed to have been redistributed uniformly after the 
stripping. The enclosed mass at r — 0.1r v i r is increased by 
22% (for / sat = 20%), 26% (for / sat = 30%), or 24% (for 
/sat = 40%), and this enhancement factor grows rapidly to- 
wards smaller radii. 

The model's features are illustrated in Figure [2] which 
shows the fractional satellite contribution to the halo mass 
profile. In each case, the fraction is fixed to / sa t at the virial 
radius and unity at r — > 0. There is a break in the curves at 
r = 0.44r v i r , which is the maximum cluster radius that can 
receive a contribution from the satellite halos. This maxi- 
mum arises from the higher virial density at z — 2 com- 
pared to the corresponding value at z = 0.3. At each radius, 
the satellite mass fraction reaches its highest value during 
stripping, when it enjoys the enhancement due to adiabatic 
compression but does not yet suffer the reduction due to 
the redistribution of baryons. At r — 0.1r v i r this fraction is 
between 60 and 90%, which makes our results at this radius 
(as seen in Figure [T] and below) relatively insensitive to the 
precise value of / sa t; higher values of / sa t push the satellite 
fraction higher towards unity, and in this limit stripping can- 
not significantly increase the density of the satellite's dark 
matter, as reflected in equation J7). 

The implications of adiabatic compression for gravita- 
tional lensing are displayed in Figure [3] which shows the 
profile of the enclosed, projected 2-D mass density. We con- 
sider the same cluster and satellite halo parameters as in the 
previous figures, and focus on the range of projected radius 
corresponding to observed cluster Einstein radii. For each 
of the observed clusters, we show the critical lensing den- 
sity versus effective Einstein radius, as a central point plus 
1-cr error ellipse. The figure is co nsistent with the factor 
of ~ 2 discrepancy highlighted by iBroadhurst &: Barkanal 
(2008) between the observed te and the median theoreti- 
cal prediction from pure dark matter simulations; here the 




Figure 1. Profiles of the enclosed 3-D mass in clusters as a func- 
tion of radius. We consider a host halo that follows an NFW 
model at z = 0.3, and assume satellites with c v j r = 4 at redshift 
Zsat = 2 that make up a total fraction / sa t of the cluster's virial 
mass. We consider a host c v j r = 6.8 and / sat = 20% (short-dashed 
curves), a host c vlr = 7.5 and / sa t = 30% (solid curves), or a host 
c vir = 8-0 an d /sat = 40% (long-dashed curves). In each case 
we show the host profile in the absence of cooling and adiabatic 
compression (bottom curve), and the final profile (top curve) of 
a host that accreted satellites that underwent gas cooling, adia- 
batic compression and stripping of their dark matter halos before 
their baryons were redistributed throughout the cluster. In each 
case / S at is normalized so that the host's dense core arises entirely 
from the satellite cores. 

typical value of r/r v jr (for the bottom curves) at the E cr ob- 
served for each cluster is smaller only by a factor of ~ 1.5, 
since our model has required us to adopt somewhat higher 
than average cluster halo concentrations. Still, these predic- 
tions in the absence of baryonic cooling for the most part lie 
well outside the 1-a error ellipses of the observed clusters. 

With the natural parameters that we have assumed for 
the satellites, adiabatic compression resolves the current dis- 
crepancy; it boosts the theoretical predictions enough to 
bring them well within the observed error ellipses. For ex- 
ample, if /s a t = 30% then for A1689, the observed (cen- 
tral) value r/r v i r = 0.072 for the effective Einstein radius 
can be compared with the predicted value (at the same pro- 
jected surface density equal to the critical lensing density) of 
r/r v - u = 0.064 (with adiabatic compression in the satellites), 
and the previous pure dark-matter prediction (i.e., without 
adiabatic compression) of r/r v i r = 0.045. For A1703, the ob- 
served r/r v i r = 0.070 can be compared with the theoretical 
r/r v - u = 0.070 (with adiabatic compression) and r/r v i r = 
0.052 (without). For C10024-17 the corresponding numbers 
are r/r v i r = 0.100 compared to r/r v i r — 0.081 (with) and 
r/r vir = 0.064 (without); for RXJ1347, r/r vir = 0.109 com- 
pared to r/r v ir = 0.100 (with) and r/r v i r = 0.086 (without). 
While the observed clusters still have slightly high Einstein 
radii compared to the typical expected cluster profile, the 
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Figure 2. Radial profile of the fractional satellite contribution 
to the enclosed halo mass. Assumptions and notations are the 
same as in Figure [l] and in particular we consider / sa t = 20% 
(short-dashed curves), 30% (solid curves), or 40% (long-dashed 
curves). In each case we show the satellite fraction in the absence 
of cooling and adiabatic compression (/P a t(r), bottom curve), the 
final fraction after stripping and baryon redistribution (based on 
equation J9}, middle curve), and the higher fraction present dur- 
ing stripping but still before baryon redistribution (equation 
top curve). 
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Figure 3. Profile of the enclosed projected surface mass den- 
sity versus projected radius. Assumptions and notations are the 
same as in Figure [l] and in particular we consider / sa t = 20% 
(short-dashed curves), 30% (solid curves), or 40% (long-dashed 
curves). For comparison, we show the observed values for four 
clusters, A1689 (open circle), A1703 (square), C10024 (triangle), 
and RXJ1347 (x). For each cluster we show its critical density 
for lensing versus effective Einstein radius, with the dot indicat- 
ing a central location and the ellipse showing the combined 1— a 
uncertainties due to the measurement errors in te and M v j r . 



theoretical scatter in c v i r together with the observational 
errors make the theoretical and observational predictions 
consistent with each other. 

As we have shown, our results depend only weakly on 
/ sa t, as long as it is within a reasonable range. The results 
also depend slightly on other assumed properties of the satel- 
lites. We illustrate this for A1689, fixing / sat = 30% and ad- 
justing the host concentration accordingly in each case. We 
find that lowering the satellite c v i r to 3 at z sa t = 2 decreases 
the predicted r/r v i r by 14%, while raising c v i r to 5 increases 
it by 12%. Assuming c v i r = 3 at z sat = 3 raises the predicted 
r/r v i r by 3%, while c v i r = 5 at z sat = 1 lowers it by 13%, 
all compared to our standard case of c v i r = 4 at z sa t = 2. 
Finally, if we assume that only 50% (rather than 100%) of 
the baryons in the galactic satellites cooled and condensed 
before their halos were stripped, i.e., in equation ([5]) we use 
half the cosmic fraction for /(,, then the predicted r/r v i r is 
reduced by 9% for A1689 and / sat = 30%. 



4 DISCUSSION 

We have demonstrated that dark matter compression due to 
baryonic cooling inside galaxy halos can in turn lead to tidal 
stripping of these galactic halos closer to the center of the 
galaxy cluster in which they reside. Even if the baryons are 
later redistributed within the cluster by feedback, a substan- 
tial effect remains due to the early adiabatic compression. 
This effect can explain the high central mass concentration 



of clusters in lensing observations. Our scenario, in which 
only the inner 10-20% of the virial radius is significantly 
modified, is consistent with weak lensing measurements at 
larger radii that find low cluster halo concentrations (e.g., 
iMandelbaum et al.ll2008T ). 

We have adopted a number of simplifying approxima- 
tions in showing the existence of the effect. Hydrodynam- 
ical simulations that avoid overcooling of the baryons at 
the cluster core are necessary in order to test our proposed 
mechanism in quantitative detail. Nevertheless, our simpli- 
fied treatment has demonstrated the general point that it is 
possible for gas physics to significantly change the Einstein 
radius of massive clusters, even without leaving a central 
baryon concentration. 

Finally, we note that ram pressure stripping of hot ac- 
creted baryons, which make up the majority (~ 80%) of the 
cluster baryons, may help to reduce the central baryon frac- 
tion. Indeed, cluster simu lations find a reduced baryon frac- 
tion at r/r v ir < 0.2 (e.g., fKravtsov et al. I l2005l : iDolag" et al.l 
2009) and match be tter the baryonic frac tion inferred from 
X-ray observations ( Vi khlinin et alj|2009l ). 
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